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I. INTRODUCTION 

The principle of superposition is at the heart of quan- 
tum mechanics. Microscopic quantum superpositions 
are routinely observed in experiments, but for macro- 
scopic objects the decoherence time is almost instan- 
taneous [TJ [2|. The famous thought experiment of 
Schrodinger [3] illustrates the problems caused by the 
transition from quantum to classical regime in the de- 
scription of physical systems of increasing size. At the 
quantum-classical boundary, systems should be isolated 
long enough, so that creation of superposition of distin- 
guishable states might be possible at least on mesoscopic 
scale. If quantum mechanics and classical mechanics are 
limiting cases of the same dynamics, the theory, which de- 
scribes the relation between them, must be nonlinear and 
stochastic [4J. Nonlinearity is crucial for the creation of 
superpositions of macroscopically distinguishable states, 
the so called macroscopic quantum superposition states 
(MQSS). These states were first suggested by Yurke and 
Stoler [5 . A small number of experimental realizations 
has been reported. Friedman et al. [6 have shown that a 
SQUID can be put into a superposition of two magnetic- 
flux states: one corresponding to a few microamperes 
of current flowing clockwise, the other corresponding to 
the same amount of current flowing anticlockwise. Brune 
et al. [7J [8] produced superpositions involving microwave 
photons. MQSS with trapped ions have been created in 
Ref. [9 . Ourjoumtsev et al. [10] also demonstrated ex- 
perimental MQSS using homodyne detection and photon 
number states as resources. 

During the last decade ultracold gases became a very 
important testbed for many predictions of condensed 
matter physics [TT] . The underlying reason for this popu- 
larity is the high flexibility of the system parameters: the 
strength of interaction, the type and strength of exter- 



nal potential(s) the particles are moving in, or the phase 
difference between two (or more) condensates. This ver- 
satility allows us to check and examine, for example, 
predicted delicate phases in different dimensions, such 
as the Bardeen-Cooper-Shrieffer [T2J [13] or the Mott- 
insulator phase [14] in 3D, the Berezinksii-Kosterlitz- 
Thouless [15j [16] phase in 2D and the Tonks- Girardeau 
phase in ID [l7l[T8]. With the advent of Bose-Einstein 
condensates (BECs) a new experimental system has be- 
come available for the investigation of macroscopic quan- 
tum superposition states. Recently, several proposals for 
generation of MQSS in BECs have been reported in the 
literature. Cirac et al. [19J have provided a method which 
involves adiabatic transfer of a two-species condensate to 
the many-body ground state. Gordon and Savage [2Q| [2Tj 
use a two-component BEC with two-body interactions 
and weak Josephson coupling between the two compo- 
nents. In Ref. [22] a scheme of creating MQSS with 
scattering light from two BECs moving with opposite ve- 
locities has been proposed. MQSS of 7r-phase imprinted 
single component BEC have been predicted to exist in a 
double- well potential [23] . 

In this study we work out a quasi- ID Hamiltonian 
model which describes two-component BEC states. In 
the two-mode approximation we solve the eigenvalue 
problems and construct MQSS (cat states) for small BEC 
systems appropriate for numerical treatment. By inves- 
tigating phase decoherence and time development we es- 
tablish signatures of cat states for larger BEC systems 
which develop into mixture states. 



Section O deals with the model, in section III 



vestigate MQSS in a double-well potential. Section IV 



we m- 
is 



left for treatment of decoherence and section IVl contains 
the conclusion and a summary. 



II. THE MODEL 
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A two-component Bose-Einstein condensate can be de- 
scribed with the following second quantized Hamiltonian 
at zero-temperature: 
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H — Hi + H2 + Hi2, 



(i) 



where 



Hi = J d 3 r ^-^-^(r)A^(r) + *t(r)V-(r)*»(r) 
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<?i 2 ^ d 3 r^ t 1 (r)^(r)^i(r)^ 2 (r). (3) 



V$(r) is the external trapping potential. The coupling 
constants ga and #i 2 are the intra- and inter-species 
atom-atom interaction strengths: 
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where mi 2 denotes the reduced mass 
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and ciij is the s-wave scattering length which can be tuned 
via Feshbach resonance. ^J(r) and ^i(r) are the bosonic 
creation and annihilation operators for the two species 
satisfying the usual commutation rules and normalized 
to the number Ni of particles of species i. 

We will assume that the external potentials have the 
following form 



Vi(r) = -m^l^y 2 + z 2 ) + V^x), (6) 



where 



V l {x) = \m i u J l x (x* + Ae- B * 2 ). 



We define 
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where <j>^ and <j>^ are the mean-field ground and first 
excited states for an atom of species i, and they are nor- 
malized to one. a[°^ and are the annihilation oper- 
ators in the ground and excited state, respectively and 
satisfy the usual commutation relations. (j)f^ and (j)^ 
can be obtained from the coupled Gross-Pitaevskii (GP) 
equations [24, 25 . The confinement of the potential is 



relatively strong in the radial directions, therefore, one 
may assume that the evolution in the weak x-direction 
decouples from that of in the strong [yz]-plane. Us- 
ing these simplifying assumptions the three-dimensional 
equation can be reduced to the following coupled, time- 
independent, one-dimensional equations [26] 
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where n\ and /j,2 are the chemical potentials and 
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foil = 2an7Vi, fo 2 2 = 2a 22 7Vift/7, &12 = &21 = 2ai 2 iV~i(l + 
+ 7)j 7 = ^2,_l/^i,_l, k = mi/m 2 , Ai = ui iX /wi,±, 
A2 = ^2, x /^i,±: A = A/clq, and B = B Oq. The distances 
are measured in units ao = y/h/m~iUJi^±. The normaliza- 
tion is such that J |0i| 2 dx = 1 and J |</> 2 | 2 dx = N 2 /Ni. 
Moreover, the relation j 2 = ft must hold if both species 
experience the same HO potential. 

According to the two-mode approximation [27] we ex- 
pand the field operator as follows 



[af 0f (a) + af 0f (a?)] x^(</, 



(12) 



where Xi,± satisfies the corresponding 2D HO 
Schrodinger equation. It is convenient to introduce 
angular momentum operators 
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where Jf , Jf correspond to the correlation between the 
two wells and Jf is the particle number difference be- 
tween the two wells for the species i. Furthermore, J 2 is 
the Casimir invariant of the SU(2) algebra: 
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By substituting (12) into gives the following result 
within the improved two-mode approximation [29] 130] 
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We measure the mass in units of the unified atomic mass 
unit. 

A similar model was introduced within a simple two- 
mode approach in Ref. [51] . 



III. MACROSCOPIC QUANTUM 
SUPERPOSITON STATES 

In this section we want to determine the domain of 
ctij's where it is possible to create cat states. We solve 



the eigenvalue problem of the Hamiltonian (15]). The 
state space can be spanned in the (Ni + l)(iV2 



1) Fock 

state basis \n^N\ — nf) 1 [n^,^ — ^2)2? wnere n f 18 
the number of particles of species i in the left well. An 
arbitrary state reads as follows 

l$>= Yl £ c nf?n ^ |nf,iV 1 -nf> 1 \n^N 2 - n£) 2 , 

=0 n£=0 
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where 
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To obtain the parameters appearing in the Hamiltonian 
( 15 ) we need to determine the mean-field ground and first 

from the time-independent 



excited states 



and 



coupled Gross-Pitaevskii equations ( 10a )-( 10b). This can 



be done in case of the ground state ^°M)y using the 
imaginary time evolution method [32] combined with the 
split-step operator technique [33 j. For getting the first 
excited state we use the shooting method [54] . 



We define P as the number N cat of superpositions of 
two distinguishable states (associated with the cat states 
or MQSS) divided by the total number N to t = [N\ + 
1)(N 2 + 1) of eigenstates ordered in sequence of growing 
eigenenergies, 



N t< 



(19) 



In practice the selection of cat states proceeds via the 
search of the N tot states (obtained by diagonalisation of 



(15)) for those members which possess two expansion co- 



efficients \c n L^ n L\ greater than 0.4. 

For our numerical investigation we fix the following pa- 
rameters: mi = rri2 = 87, uji^± = U2 y ± = 2tt x 710 Hz, 
Ai = A 2 = 0.2, A = 50, B = 10. This means 
that our exploratory investigation for two-component cat 
states will be carried out for the case of a quasi one- 
dimensional (cigar-like) 87 Rb Bose gas where the two 
components correspond to two hyperfine states of the 
rubidium atom. Such two-component BEC systems are 
routinely produced in laboratories [35] with particle num- 
bers Ni ~ N2 ~ 10000. But simulating such a great 
number of atoms in the gas is beyond the capacity of 
today's computers. We think, however, that features of 
cat states to be presented here for lower particle num- 
bers Ni and N 2 will be valid also for larger Ni values. 
Therefore we shall carry out the calculations at particle 
numbers N\ , N2 ~ 1 — 20 and vary the scattering lengths 
in the possible domain which can be achieved by using 
the Feshbach resonance method [56] . 

In Fig. [T] the probability P of cat state formation is 
shown as function of an and a\2 with fixed parameters 
<222 = 2 nm, and particle numbers N\ = 4 and AT 2 = 30. 
We see an almost symmetrical triangle which is cut at 
a\2 = where there is no two-component cat states. 
This is reasonable because the components are decou- 
pled and, at the most, one-component cat states can be 
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FIG. 1: (Color online) P as function of an and an in nm; 
fixed parameters: Ni — 4, N 2 — 30, a 2 2 — 2 nm. 



formed with such parameter values. We observed how- 
ever also that there is no two-component cat states when 
the sign of CL22 is changed to negative (at the above fixed 
parameters). 

Therefore, it is reasonable to study the formation of 
cat states as function of the particle numbers N\, N2. 
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FIG. 2: (Color online) P as function of N\ and N 2 ; fixed 
parameters: an — 2 nm, a± 2 — 4 nm, a22 — —2 nm. 

In Fig. [2]P is shown as function of N\ and N 2 at posi- 
tive an,ai2 and negative CL22 values. A sharp horizontal 
border can be observed from N2 > 10, independent of 
N±. This result explains the former finding of no cat 
state when a 2 2 has been changed to negative values. We 
numerically observed also that this horizontal phase sep- 
aration border (between cat states and no cat states) can 
be raised to larger N2 values by raising also the value of 



In Fig. 3(a) all the scattering lengths have been fixed at 
positive values. It can be seen that the probability P of 
forming cat states is lowered at growing particle numbers 
(decline is slower for component with smaller scattering 
length). But this means only that greater particle num- 
bers involve more N tot states and thus a smaller fraction 
becomes cat states. This can be visualized if the result 
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FIG. 3: (Color online) Fixed parameters: an = 2 nm, ai2 — 
4 nm, a22 = 3 nm. (a) P as function of Ni and N 2 . (b) N cat 
as function of Ni and N 2 . 



of the same calculation is shown not for the probability 
P but for the number N cat of cat states. In Fig. 3(b) 



we 



see that the number of cat states are growing up to 140 
at particle numbers N± = N2 = 20 which means however 
total number of states N tot = 441, and therefore gives for 
P only the value 0.3 which has been represented by blue 



color in Fig. 3(a) _ 

In Fig. |4(a)| the intra-species scattering lengths an, 022 
are positive while the inter-species scattering length a\2 
is negative. According to the color bars we se e tha t P 
decreases here more strongly compared to Fig. |3(a)| but 
the nu mber o f cat states is also decreasing, as can be seen 
in Fig. |4(b)| where N cat is visualized for the same set of 
parameters. 

In Fig. |5(a)| and Fig. 5(b) all scattering lengths take 
negative values. We see that for larger particle numbers 
there are no cat states at all. 



IV. DECOHERENCE 

We study the effect of dissipation. The dominant 
source is the thermal cloud of noncondensed atoms. 
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FIG. 4: (Color online) Fixed parameters: an = 2 nm, ai2 — 
—4 nm, a22 = 3 nm. (a) P as function of Ni and iV 2 . (b) 
N cat as function of Ni and AT 2 . 



FIG. 5: (Color online) Fixed parameters: an = —2 nm, ai2 = 
—4 nm, a 2 2 = —3 nm. (a) P as function of N± and iV 2 . (b) 
N cat as function of iVi and N2. 



Other possible sources are ambient magnetic fields, dif- 
ferent scattering lengths, three-body losses [37 J. We will 
only investigate the effect of elastic two-body interactions 
between the noncondensed and condensed atoms. This 
type of interaction preserves the number of condensed 
atoms and leads to phase-damping which dominates for 
MQSS HQ. 

The most general quantum mechanical Markovian 
master equation has the form [38] 



dt 



(20) 



where p is the density matrix of system states, £ a 's are 
the Lindblad operators containing system-environmental 
interaction, and {•, •} is the anticommutator. The Lind- 
blad formalism is valid when the coupling between the 
system and environment is weak [39] . We use the follow- 
ing Lindblad generators to describe the phase-damping: 



Ci = 

c 2 = 

U = 



/vmf, 

K2U2 , 



(21a) 
(21b) 
(21c) 
(21d) 



where K\ and K2 are the phase-damping rates, which 
are proportional to the number density of noncondensed 
atoms. 

The dynamics of the master equation is studied in 
terms of stochastic trajectories of state vectors [40| |4T]. 
In this approximation the density matrix of the system 
is written as the average of the stochastic wavefunctions: 



MC 



rt*) = McEI $(i) (*))( $(i) (*) 



(22) 



where |$W i s the zth stochastic trajectory of the state 
vector, and MC is the number of the simulations. In 
the evolution of the system two distinct elements can 
eventuate for small At. Either one of the quantum jumps 
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occurs (with the probablity dp a = (& a C a ) At) which can 
be represented as an action of the corresponding Lindblad 
operator on the state of the system 



(23) 



or the wavefunction is evolved by a non-Hermitian ef- 
fective Hamiltonian (with the probablity 1 — ^2 a dp a ), 
which has the form: 



H eff H 



ik 



(24) 



k=l 



The trajectories must be normalized after every time 
step. The Monte-Carlo approach has many advantages. 
Calculating a set of trajectories is more economical, than 
solving the master equation with a Runge-Kutta method. 
Moreover, a single-stochastic trajectory can be inter- 
preted as a single realization of an experiment manip- 
ulating a unique quantum system [42] . 

Our calculation will be performed with parameter set 
an — 2 nm, a 12 = 4 nm, (122 = 3nm, Ni = 5, 
N2 — 20 and phase damping rates K\ = K2 = 0.005. 
The density matrix is obtained by averaging MC = 
5000 trajectories. Initially, among the 126 e igenstates 
there are about 60 cat states as Fig. |3(b)| suggests. 
From these we shall investigate only two which can be 
written to a good approximation as (|1, 4) 1 |20, 0) 2 + 
14,1)^0,20)2)/^ (called the 61th excited state) and 
(12,3)^7,13)2 + 13,2)^13,7)2)/^ (the 63th excited 
state) using the standard notation of the expansion (17). 

Fig.[6]shows the magnitude of the density matrix of the 
61th excited (cat) state at different times. At initial time 
t = there are two off-diagonal elements which are equal 
to the two diagonal ones. While the diagonal elements do 
not change significantly in time, the off-diagonal ones fall 
off gradually, and the state becomes a classical mixture 
after 2 second. A similar behavior (not shown) can be 
observed for the 63th cat state (and for all the other cat 
states). 

To study finer details of the decoherence effect of en- 
vironment on cat states we have calculated both Tip 2 
and the magnitude C of the average of the two maxi- 
mal off-diagonal elements as function of time for the se- 
lected 61th and 63th excited states. These two quantities 
are characteristic of time development of the decoherence 
process, measuring formation of statistical mixture from 
the initital pure cat state. 

The results are shown in Figs. [7] and [8] and we can ob- 
serve that the decoherence is more effective for the state 
where the difference AA^ = \N^ — Nf-\ of particle num- 
bers of the two components i = 1, 2 in the left and right 
potential well is greater. For example, in case of the 
61th excited cat state, ANi = 3 and AA2 = 20 which 
is greater than the corresponding values AAi = 1 and 
AiV 2 = 6 of the 63th state, and thus Trp 2 of the 61th 
state approaches more quickly the value 0.5, characteris- 
tic of the statistical mixture than that of the 63th state. 
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FIG. 6: The magnitude of density- matrix at different times 
(61th excited state): (a) initial state, (b) t = 0.5 sec, (c) 
t = 2 sec. Fixed parameters: an — 2 nm, a±2 — 4 nm, 
a 22 = 3nm, Ni = 5, A 2 = 20, m = k 2 = 0.005. 



A similar behavior can be observed when the time devel- 
opment of the average C of the two maximal off-diagonal 
density matrix elements are studied in Figs. [7] and [§J 
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FIG. 7: The 61th excited state: (|1, 4> x |20, 0> 2 + 
|4, 1) 1 |0, 20) 2 )/v / 2- Fixed parameters: an = 2 nm, an — 
4 nm, a 22 3nm, iVi = 5, A 2 = 20, «i = k 2 = 0.005. 
Time is measured in units of second, (a) Trp 2 as function of 
time, (b) The magnitude C of the average of the off-diagonal 
elements as function of time. 
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FIG. 8: The 63th excited state: (|2, 3) 1 |7, 13> 2 + 
|3, 2) 1 1 13, 7) 2 )/v / 2- Fixed parameters: an = 2 nm, ai 2 — 
4 nm, a 22 3nm, iVi = 5, A 2 = 20, «i = k 2 = 0.005. 
Time is measured in units of second, (a) Trp 2 as function of 
time, (b) The magnitude C of the average of the off-diagonal 
elements as function of time. 



V. CONCLUSION 

We have studied macroscopic quantum superposition 
(cat) states formed in two-component Bose-Einstein con- 
densates described within the two-mode approximation. 
An analysis of the limitations of this approximation can 
be found in Ref. [43 . However, the two- mode approx- 
imation has been ascertained to describe many aspects 
of double well experiments [44]. The model has been 
applied to BECs where the two components correspond 
to two internal hyperfme state of 8T Rb atoms. Because 
of limited computational capacity, relatively modest par- 
ticle numbers between 1-20 has been investigated, but a 
systematic check revealed that the main features of form- 
ing the cat states depend mostly on the actual interaction 
strengths (characterized by scattering lengths a^) be- 
tween atoms, and not on the actual numbers A^, i = 1, 2. 
Therefore we have performed exploratory calculations in 
wide range of values of scattering lengths a^-,z, j — 1,2, 
both in the attraction and repulsion domain, and pre- 
sented the results in case of some selected examples ex- 
hibiting the most important features of formation, dura- 
tion, and ceasing of MQSS or cat states. 



We have found that, in general, positive inter-species 
scattering lengths a^,i = 1,2 enhance the probabil- 
ity P = N ca t /N tot of two-component cat state forma- 
tion. Small or zero intra-species scattering lengths ai2, 
of course, lessen this probability, as the system of two- 
component BEC actually decouples into two independent 
one-component BECs where, however, formation of one- 
component cat states is possible. In case of positive inter- 
species scattering length a\2 also the individual number 
N cat of two- comp onent cat states is growing as it can be 
seen in Fig. 3(b)| b ut ne gative a\2 hinders the formation 
of cat state as Fig. |4(b)] shows. Finally, in the case when 
all the scattering lengths take negative values, there is 
only a limited range of two-component BEC cat state 
formation, studied both in terms of probability P or the 
individual cat state number N cat , as Figs. 5(a)| and |5(b)| 
exhibit. 

Effect of environment has been modeled by phase- 
damping which is the result of elastic interaction between 
condensed and non-condensed atoms and preserves the 
particle numbers 7V l5 N2. Using the corresponding Lind- 
blad master equation (20) we have found that the most 



important parameter which gives rise to the decoher- 
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ence effect of such type is the particle number difference Acknowledgment 
ANi = \N t L - N t R \,i = 1,2 between the left and right 
potential well. As Figs. [7] and [8] suggest the decoherence 
occurs more quickly for two-component cat states where 
the number difference ANi is greater. 

The influence of asymmetric wells, and other effects 
of the environment on the formation of cat states needs The authors thank Dr. Daniel Schumayer for reading 
further research. the manuscript and providing valuable comments. 
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